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^ (57) Abstract: Techniques are described for estimating the propagation velocity through a surface acoustic wave sensor. In partic- 
ular, techniques which measure and exploit a proper segment of phase frequency response of the surface acoustic wave sensor are 
described for use as a basis of bacterial detection by the sensor. As described, use of velocity estimation based on a proper segment 
of phase frequency response has advantages over conventional techniques that use phase shift as the basis for detection. 
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ESTIMATING PROPAGATION VELOCITY THROUGH 
A SURFACE ACOUSTIC WAVE SENSOR 

CROSS-REFERENCE TO RELATED APPLICATIONS 

5 This application claims priority to U.S. Provisional Patent Application 

No. 60/533,177, filed December 30, 2003. 

« 

TECHNICAL FIELD 

The invention relates to surface acoustic wave (SAW) sensors and, more 
10 particularly, techniques for analyzing and interpreting the output of a SAW sensor. 

BACKGROUND 

* 

Chemical and biological testing is commonly used to test for the presence or 
absence of chemical or biological agents. Testing for the presence of chemical or 

15 biological agents in blood, food or other materials is often performed to ensure safety 
or to facilitate diagnosis of medical conditions. For example, testing is used to identify 
chemicals, bacteria or other agents in blood samples taken from medical patients, 
laboratory samples developed for experimental purposes, food samples, or the like. In 
addition, chemical and biological testing is used to test for medical conditions such as 

20 pregnancy, diabetes, bacterial infection, and a wide variety of other conditions that may 
affect the patient's chemistry or biology. 

One type of sensor that has been developed for chemical or biological sensing 
capabilities is a surface acoustic wave (SAW) sensor. One example of a SAW sensor is 
a Love mode shear-horizontal surface acoustic wave (SH-SAW) sensor. A SH-SAW 

25 sensor includes four main components: 1) a piezoelectric substrate; 2) an input inter- 
digitated transducer (IDT) on the substrate, which is used to excite an acoustic wave 
based on the piezoelectric effect; 3) an output IDT on the substrate, which receives the 
transmitted acoustic wave and generates electrical output by exploiting the piezoelectric 
effect; and 4) a wave-guide layer over the IDT's, which converts SH-type waves into 

30 waveguide Love modes for transmission from the input DDT to the output IDT. The 
presence of one or more materials on the surface of the SH-SAW sensor affects wave 
propagation through the waveguide layer in response to the presence of bacteria or 
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other agents on the surface of the sensor, which facilitates detection of bacteria or other 
agents. 

SUMMARY 

i 1 

% 

5 In general, techniques are described for estimating the propagation velocity, or 

equivalently, for estimating a time delay through a surface acoustic wave sensor. In 
particular, techniques which measure and exploit a proper segment of phase frequency 
response of the surface acoustic wave sensor are described for use as a basis of 
bacterial detection by the sensor. 

* 

10 In one embodiment, the invention provides a method comprising identifying a 

segment of phase frequency response of a surface acoustic wave sensor, and estimating 
a time delay associated with wave propagation through the surface acoustic wave 
sensor based on the identified frequency response. 

In another embodiment, the invention provides a computer-readable comprising 

1 5 instructions that when executed in a processor identify a segment of phase frequency 
response of a surface acoustic wave sensor, and estimate a time delay associated with 
wave propagation through the surface acoustic wave sensor based on the identified 
frequency response. 

In another embodiment, the invention provides a system comprising a surface 

20 acoustic wave sensor, a sensor analyzer to receive output of the surface acoustic wave 
sensor, and a processor to receive input from the sensor analyzer, identify a segment of 
phase frequency response of a surface acoustic wave sensor, and estimate a time delay 
associated with wave propagation through the surface acoustic wave sensor based on 
the identified segment of phase frequency response. 

25 The invention may be capable of providing one or more advantages. In 

particular, use of change of propagation velocity of a surface acoustic wave sensor, as 
described herein, can improve detection of bacteria via the sensor, relative to 
conventional techniques that use phase shift as the basis for detection. Moreover, use 
of an estimated propagation velocity of a surface acoustic wave sensor as the basis for 

30 detection may allow for detection of bacterial concentration. 

The details of one or more embodiments of the invention are set forth in the 
accompanying drawings and the description below. Other features, objects, and 
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* 

advantages of the invention will be apparent from the description and drawings, and 
from the claims. 

BRIEF DESCRIPTION OF DRAWINGS 

5 FIG. 1 is a perspective view illustrating an exemplary surface acoustic wave 

■ 

(SAW) sensor that may be used in one or more embodiments of the invention. 

■ 

FIG. 2 is a block diagram illustrating a system according to an embodiment of 
the invention. 

FIG 3 is a flow diagram illustrating a technique according to an embodiment of 

10 the invention. 

FIG. 4 is a graph illustrating an example phase frequency response of a SAW 

■ 

sensor. 

FIG. 5 is a depiction of an exemplary Lab View screen associated with 
execution of the techniques described herein. 
1 5 FIG. 6-3 1 are various graphs illustrating techniques that can be used in 

accordance with the invention, and various desirable characteristics that can be 
observed relative to the characteristics observed by conventional techniques. 

DETAILED DESCRIPTION 

20 FIG. 1 is a perspective view illustrating an exemplary SAW sensor 10 that may 

be used in accordance with an embodiment of the invention. SAW sensor 10 may 
comprise any of a wide variety of SAW sensors. SH-SAW sensors are typically 
constructed from a piezoelectric material with a crystal-cut and orientation that allows 

* 

the wave propagation to be rotated to a shear horizontal mode, i.e., parallel to the plane 
25 defined by the waveguide, resulting in reduced acoustic damping loss to a liquid in 
contact with the detection surface. Shear horizontal acoustic waves may include, e.g., 
thickness shear modes (TSM), acoustic plate modes (APM), surface skimming bulk 
waves (SSBW), Love-waves, leaky acoustic waves (LSAW), and Bleustein-Gulyaev 
(BG) waves. 

30 In particular, Love wave sensors may include a substrate supporting a SH wave 

mode such as SSBW of ST quartz or the leaky wave of 36°YXLiTa0 3 . These modes 
may preferably be converted into a Love-wave mode by application of thin acoustic 
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guiding layer or waveguide. These waves are frequency dependent and can be 
generated if the shear wave velocity of the waveguide layer is lower than that of the 

piezoelectric substrate. 

In one example, sensor 10 comprises a Love mode shear-horizontal surface 

5 acoustic wave (SH-SAW) sensor. 

SAW sensor 10 includes a substrate 12 which typically comprises a 

piezoelectric material. SAW sensor 10 also includes an input inter-digital transducer 

< 

(IDT) 14 on substrate 12, which is used to excite an acoustic wave based on the 
piezoelectric effect. In addition, SAW sensor 10 includes an output IDT 16 on 

1 0 substrate 1 2, which receives the transmitted acoustic wave and generates electrical 

output by exploiting the piezoelectric effect. A wave-guide layer 18 is formed over the 
IDT's 14, 16. Wave-guide layer 18 converts SH-type waves into waveguide Love 
modes for transmission from input IDT 14 to output IDT 16. 

A layer of material, such as a layer of antibodies is applied over waveguide 

15 layer 18. In operation, a fluid being tested for the presence of bacteria is brought into 
contact with waveguide layer 1 8. If bacteria is present in the fluid, the bacteria attaches 
to the antibodies on waveguide layer 18 and thereby affect wave propagation through 
waveguide layer 18. Accordingly, analysis of the wave propagation through waveguide 
layer 18 of SAW sensor 10 allows for bacterial detection or detection of other agents 

20 that may interact with a material coated on waveguide layer 1 8. 

In some cases, SAW sensor 10 includes a plurality of sets of injput and output 
IDT's. For example, SAW sensor 10 may include a first input IDT 14 and a second 
input IDT 15 which respectively correspond to first output IDT 16 and second output 
IDT 17. In that case, first input and output IDT's 14, 16 comprise the active portion of 

25 SAW sensor 10, and second input and output IDT's 15, 17 comprise the reference 

portion of SAW sensor 10. Different types of antibodies may be coated on the surface 
of waveguide layer 18 between first input and output IDT's 14, 16 and second input and 
output IDT's 15,17, such that only the bacteria of interest bonds to the antibodies 
between the active portion corresponding to first input and output IDT's 14, 16. In that 

30 case, the reference portion corresponding to second input and output IDT's 15, 17 

allows for reference measurements which can account for temperature variance affects, 
or the like, which could otherwise affect wave propagation through SAW sensor 10. 
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SAW sensors are commonly used for bacterial detection but may be designed 
for detection of any of a wide variety of other chemical or biological agents. 
Accordingly, different materials maybe coated on the waveguide layer of SAW sensor 
10 in order to facilitate detection of various chemical or biological agents. Waveguide 

5 materials may preferably be materials that exhibit one or more of the following 

properties: low acoustic losses, low electrical conductivity, robustness and stability in 
water and aqueous solutions, relatively low acoustic velocities, hydrophobicity, higher 
molecular weights, highly cross-linked, etc. In one example, SiC>2 has been used as an 
acoustic waveguide layer on a quartz substrate. Examples of other thermoplastic and 

10 crosslinked polymeric waveguide materials include, e.g., epoxy, 

polymethylmethacrylate, phenolic resin (e.g., NOVALAC), polyimide, polystyrene, 

etc. 

* 

In general, the presence of a particular material on the surface of the SAW 
sensor affects wave propagation through the waveguide layer in response to the 
1 5 presence or absence of another material passing over the waveguide layer, which 

facilitates detection of the material passing over the waveguide layer. Accordingly, 
materials coated on the waveguide layer may be selected to attract, trap, bond with or 
otherwise attach to materials suspended in a fluid that flows across waveguide layer 18. 
In this manner, SAW sensor 10 facilitates detection of bacteria or other biological or 

20 chemical agents. 

In order to generate Love mode surface wave, the shear wave velocity of 
waveguide layer 1 8 is typically lower than that of substrate 12. In that case, the 
acoustic energy will be trapped to near the sensing surface of SAW sensor 10. Love 
mode SH-SAW sensors, in particular, generally have high sensitivity to surface 

25 perturbations. Mass loading can change the surface conditions of a Love mode SH- 
SAW sensor. Therefore, measuring the change of the propagation velocity or resonant 
frequency can be used to quantitatively detect mass loading, and thereby detect the 
presence of a chemical or biological agent. 

One of the most important performance indices for an SH-SAW sensor used as 

30 a detector is the mass sensitivity, i.e., how sensitive is the sensor to the loading mass on 
its surface. For analyzing the mass sensitivity, two indices have been established: 
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S J,=-1 lim -*L and S f m =± lim Af 



Vq Am g ->0 Am 3 f 0 Am 8 -+0Am s 



where * = and f = y ™°" a ™ dMss™* . Vo , f <°> onance and v,, f < e 1 > onanc9 are 

■ 

the propagation velocities and resonant frequencies of the sensor without and with 
5 surface perturbation resulting from an infinitesimal thin rigid loading layer with mass 

■ 

Am s =pme» where p m and 6 are the density and thickness of the loading layer. 

Sj{, « — 5 m , where v g is the group velocity, and for Love mode SH-SAW 

v ° ^— > 1 , where the notation v 0 ((o) as well as v g (co) is used to emphasize that the 
v g (a>) 

velocities are frequency dependent in the dispersive cases. This means that for Love 

1 0 mode SH-SAW sensor, using — for detecting the anomalies of the boundary 

v o 

conditions of the surface of the sensor should be better than using ^ . Accordingly, in 

fo 

accordance with the invention, —may be used as a basis of a detection indicator for 
SAW-sensor 10, or the like. 

« 

Conventionally, for a Love mode SH-SAW sensor with a triple transit echo 

> 

1 5 (TTE), the phase frequency response was given by: 



. x 1 + p Leo 

(01©) = , 

Vl ' 1-pv(co)' 



or equivalently: 



^ ' 1-p v(f) ' 



where f and a> represent frequency and angular frequency, L is the distance between the 
20 center of the input and center of the output IDT, p=a and a is the reflection coefficient 
of the input and output DDTs. Since P is small and can be determined by frequency 
response, without loss of generality, p is assumed to be zero. Thus, for phase frequency 
response, one only needs to be concerned with the case that does not include TTE. The 
phase frequency response is therefore: 
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However, e^<" = e^'* = e-^WO+k) , where k=0, *1, ±2, +3 etc. Thus, the 
phase frequency response <p(f) is multi-valued and cannot be uniquely determined by 
the value of o Mf) . This problem is called phase ambiguity. Only a so-called "main 
5 value" of (p(f) can be determined from e Wf > . Since 

exp{-j2nft(f )} = exp{-j2K(fi(f ) - [ft(f )])} , 

» 

where [x] is the integer part of x, i.e. [x] is the maximum of the integers, which are less 
or equal to x, the main value of the phase response becomes; 

♦(f) = -27i(fi(f ) - [ft(f )]) + 7i (in radian) or ftf ) = -360(ft(f ) - [fr(f )]) + 1 80 (in 

1 0 degree). 

Thus, -7r<<|>(f)£7r (in radian) or - 1 80<<j>(f)< 1 80 (in degree). 

For (p(f) = -2nf-, 



15 



v 



dv lv=V o = 27lf ^' 



0 

t 

Therefore, an approximation of the phase change is given by: 



A<p = 2nf — ^-Av. 

v o 



Therefore, 



27if-t-Av 

A(p = v p _ Av 

<Po " _ 2jlf JL " v 0 



Unfortunately, however, this derivation does not hold for 

20 <M f) = - 2jt (ft- [ ft )) + n . If [f-Ll = [f —L-], then |i = -^. Otherwise, 

v v v 0 v 0 + Av $0 v 0 

_L A v + [f— J^_]_ rf _L, 
A4> v o Vq+Av v 0 ^ 

*o __L +[f _Lj ""v7 - 

v o v 0 
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The condition [f— ] = lf — ^— ] implies f 

v 0 v 0 +Av 



L L 

v 0 v 0 +Av 



<1. Let f 



L 



v 0 v 0 + Av 



= 1, 



i.e. f 



LAv 



Vo(v 0 + Av) 



= 1 . For Av < 0 , it follows that: 



v 0 (v 0 +Av) fL fL + v 0 

2 

This means that when lAvl ^ - V ° , [f — ] * [f — ^-r] . For example, If 

fL + v 0 v 0 v 0 +Av 



5 f^l03MH,v 0 =4000 m/s and L=8.8mm, then [f — ] * [f — for |Av| ;> 17.5747 m/s. 

v 0 v 0 + Av 11 

Therefore, when [f — 1 * [f — - — ] , A4> or — may not be a proper indicator function. 

v 0 v 0 +Av b 0 

For this reason, the invention provides techniques for estimating the 
propagation velocity through a SAW sensor based on an identified segment of phase 
frequency response. The propagation velocity is simply related to the time delay 
10 : through the SAW sensor, as described below. In accordance with the invention, a 
method may include identifying a proper segment of phase frequency response of a 
surface acoustic wave sensor, and estimating a time delay associated with wave 
propagation through the surface acoustic wave sensor based on the identified frequency 
response. 

15 FIG. 2 is a block diagram illustrating a system 20 that includes a SAW sensor 

10, and a sensor analyzer 21 that obtains measurements from SAW sensor 10. System 
20 also includes a processor 22 that interprets the output of SAW sensor 10. In other 
words, sensor analyzer 21 receives output from SAW sensor 10 and provides input to 
processor 22 so that the output of SAW sensor 10 can be interpreted. 

20 Processor 22 receives input from sensor analyzer 21, which comprises 

measurements associated with wave propagation through SAW sensor 10. Processor 
22 then determines whether SAW sensor 10 has detected the presence of particular 
bacteria or other material for which SAW sensor 10 is designed to detect. Processor 22 
executes instructions to perform various techniques and functions ascribed to the 

25 processor herein. Although the invention is not limited in this respect, SAW sensor 10 
may be housed in a cartridge, or the like, and may be electrically coupled to sensor 
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analyzer 21 via insertion of the cartridge into a slot. Processor 22 may be housed in the 
same unit as sensor analyzer 21 or may be part of a separate unit or separate computer. 

Processor 22 may also be coupled to memory 24, which stores a propagation 
velocity routine 25 consistent with the teaching of this disclosure. Alternatively, 
5 propagation velocity routine 25 may be implemented by hardware within processor 22. 
In any case, processor 22 executes propagation velocity routine 25 in order to estimate 
a time delay associated with wave propagating through the surface acoustic wave 
sensor based on the identified frequency response, as described herein. 

By way of example, processor 22 may comprise a general-purpose 

1 0 microprocessor that executes software stored in memory 24. In that case, processor 22 
may be internally housed in a specifically designed computer, a general purpose 
personal computer, workstation, handheld computer, laptop computer, or the like. 
Alternatively, processor 22 may comprise an application specific integrated circuit 
(ASIC) or other specifically designed processor. In any case, processor 22 executes 

1 5 propagation velocity routine 25 in order to estimate a time delay associated with wave 
propagation through the surface acoustic wave sensor 1 0 based on the identified 
frequency response, as described herein. 

Memory 24 is one example of a computer readable medium that stores 
processor executable software instructions applied by processor 22. By way of 

20 example, memory 24 may comprise random access memory (RAM), read-only memory 
(ROM), non-volatile random access memory (NVRAM), electrically erasable 
programmable read-only memory (EEPROM), flash memory, or the like. Propagation 
velocity routine 25 such as one of those mathematically described below, are stored in 
memory 24 and may form part of a larger software program used for analysis of the 

25 output of SAW sensor 10. For example, propagation velocity routine 25 may be a sub- 
routine programmed within a LabView software platform, which is described in greater 
detail below. 

FIG 3 is a flow diagram illustrating a technique according to an embodiment of 
the invention. As shown in FIG .3, fluid is brought into contact with a surface of SAW 
30 sensor 10 (31). The fluid may include samples of material for which bacterial testing is 
needed. For example, the fluid may be allowed to pass over the surface of waveguide 
layer 18 of SAW sensor 10 (FIG 1), which includes antibodies that react to the bacteria 
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of interest. SAW sensor 1 0, for example, may be housed in a cartridge that defines a 
fluid path over the surface of waveguide layer 18 such that a fluid can be introduced 
into the cartridge and allowed to pass over waveguide layer 1 8- via the fluid path. 

Processor 22 receives measurements taken from SAW sensor 1 0 by sensor 
5 analyzer 21, and applies propagation velocity routine 25 stored in memory 24. In doing 
so, processor 22 identifies a segment of phase frequency response of SAW sensor 10 
(32), and estimates a time delay associated with wave propagation through SAW sensor 
10 based on the identified frequency response (33). In some embodiments, an 
estimated propagation velocity derived from the estimated time delay can be used to 
1 0 identify the concentration of the bacteria or other material in the fluid (34). 

Numerous mathematical techniques for estimating the time delay associated 
with wave propagation through the surface acoustic wave sensor will now be discussed. 

t = - refers to the time delay between the center of the input and the center of the 
v 

output IDT separated by distance L. In the dispersive case, the time delay is also a 

■ 

1 5 function of the frequency. That is, t(f ) = ^ . It follows from the definition of <|>(f) (in 



degree) that: 



Mf) = [fT(f>]-^ + 0.5. 



For any given running frequency fo (also referred to as the operating 
frequency), there are two phase inflection frequencies fi and f 2 proximate to a running 
20 frequency f 0 such that: 

(1) fjSfo <f 2 

i 

(2) 4ft ) = 180 , ((>(f 2 - 0) = -180 and <t>(f 2 ) = + 0) = 180 , 

(3) fitft) = [fix(fi)J, f 2 *(f 2 ) = l»2t(f 2 )] and [f 2 x(f 2 )] = [Wf|)]+1, and 

(4) For any f, £ f < f 2 ,[ft(f)] = fttft)] , 

25 FIG. 4 is a graph illustrating an example frequency response of a SAW sensor. 

Labeled in FIG. 4 is an exemplary a running frequency f 0 , a typical first phase 
inflection frequency fj and a typical second phase inflection frequency f 2 . Again, for 
any f, s f < f 2 , [ft(f)] = IVft)] . Therefore: 
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Hf) = 1^)1-^ + 0.5, fbrf 1 sf<ff a . 
Taking the derivative with respect to f for both sides of the above equality, results in: 



dx(f) 



df 1 ' 360 df 

That is, 

+ 

dt(f ) . 



t (f) = -J_iie),forf 1 v £ f<f 2 . 



1 ... 1 1d(|)(f) - . . . 



df 

> 

which is one basic differential equation for calculating the time delay x(f). 

If the initial condition is: i(f) | f=foo = x(f 00 ) , where f 00 € [f 1f f 2 ) and arbitrary, then 

the solution of the basic differential equation for calculating the time delay t(f) 
becomes: 

io T(f ) = Y, (foo) __l_l w+ ^_i Woo) . 

The solution can be easily derived from: 



and 



foot(f 0 o) = [M(fi)l-^ + 0.5. 



1 5 The basic differential equation may be used to prove the following properties of x(f). 

Proposition 1 : If for f e [fn,f 22 ) Q Vvh) > MO is linear wit h respect to f, then 
x(f) is constant and 

T(f) = -3i0 ^-fi ( i fl1> ' 1 6 t,11 ' ,22) S Pl ' f,) ' 
Proof of proposition 1 : 

20 Suppose f 00 e [fn,f 22 ) ■ *(0 = --i-^) is the solution of the basic differential 

360 

equation: 

f^r + x(f) = -;^r^Sr > in I'll. '22) because the left side 

df 360 df 
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f M) +T(f) = t( f ) = __l_^ oo) and the right side -£<jB---|-4cW due 

* 

to the linearity of <Kf) in [fn,f22) • It is also from the linearity of <f>(f) in[fn,f22) that 
Woo/ = — ; * • 

f 22~ f 11 

Proposition 2: Suppose <j>(f) is differentiate in (f-i f f 2 ) and right differentiate 
5 for f| . Then: 

i(n o — L is the first order estimation of t(f). 
w 360 df 

Proof or proposition 2: 

■ 

In the neighborhood of f 0 o e [f 1f f 2 ) , 

*(0«*(foo) + Woo)(^^oo). 
10 Therefore, in this neighborhood <|>(f) is approximately linear with respect to f and 

W) « Woo) • From the proposition l,x(f) « -^^*r ' That is ' * (f) = " 3io ^dP iS ^ 

first order of estimation of t(f). 

Algorithms or routines for estimating the Time Delay t(f) will now be 
discussed. From proposition 2, the following algorithms can be used as part of 
1 5 propagation velocity routine 25 to obtain the estimate of t(f). 

(0001] Algorithm 1 : f 00 - U and %(U ) = {ft ) « 1 <M f 2 ) - W1 ) ^ 1 



360 1 360 f 2 - fi f 2 - ft 

Using the equation: 



T<f) = ^x(foo)-^>)+^Woo>. 



one obtains: 



20 t(f)= ii_l__JLM + M. 

11 f f 2 360 f f 

« 

In particular: 

-„ ,_Ji__J 1 Wo) , 0.5 

(0, "'o f 2 -fi"360 f 0 + f 0 ' 

where fo is running frequency and is also denoted by fn, n in this description. 
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Algorithm 2: f 0 o = » e.g. <fr(f*) = ±90 . Around f# take 10 frequencies 
f<V < 2 >,....f < 10 > and a linear regression for <f (1 \«Mf< 1 >)),(f< 2 >,)Mf< 2 >)) (f (10, .«f (10> )). 

1 

The slope of this line is denoted by $(f . ) and x(f * ) = -— - <KM . .Therefore, 

360 

t(f) = — <ji(f . ) — — - flf ) + — - *(f. ) 

360 f 360 f w ' 360 f 

5 Algorithm 3 : f 00 = f 0 - . Around f 0 take 1 0 frequencies f (1) ,f w . . f (10) 

and a linear regression for <f (1) ,<Kf (1) )Mf (2 \4'(f (2> )) (f (10, ,<Mf <10) )). The slope of this 

line is 

* 

1 

$(f 0 ) and the first order estimation of t(f 0 ) is T(f 0 ) = — — ■ Wo ) • 

360 

Algorithm 4: Since 

10 i(0-Ywa > )-i>+i|4fa)--i^4cr»)-^Jw + i>«) 
holds for any f 00 e ft , f 2 ) . Regarding foo as a variable and taking integral with respect 

to foo for both sides of the above approximate equation, provides 

f fx(f)dfoo m --l-l f ff 00 Wf 00 )df 00 --i- f /l W f)df 0 o + -j- ) 7 «f 00 )<%> , 
f 360 f f t 360 ^ f 360 ^ f 

which is equivalent to 

15 x(f)(f 2 -f 1 )«-^l(f 2 *(f 2 )-f^(f 1 )-|[^f 00 )df 00 )-^l4>(f)(f 2 -f l)+ -l_|l W f 0() )df 00 



i.e., 

t(f) . _ jLitmizhm _ jLi W)+ _2_i_L_'f Wfoo)d , 00 

1 ' 360 f f 2 -U 360 f * ' 360 f f 2 -f, ^ °°' 00 

Because Wi ) = 180 , ()>(f 2 - 0) = -180 , 

21 fj-f, 360 f W ' 360 f f 2 -fi ^ 001 00 



20 



1 fi 



tttW + ^ + -~t . 1 , /«Mf<)0 ) df oo • 



ff 2 -fi 360 f" f isof^-fw, 

In particular, 
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T(fo)= ^^ _ 3i?fo- Wo,+ v + ^^^r ,oo)d,o °- 

Comparing the i(f) in algorithm 1 with that of algorithm 4, shows that algorithm 

« 

4 adds one more term ^-J— v ] Woo )<*foo • If the phase function <K0 is symmetric 

180 f f 2 - Ti fl 

about the point ((fi+f 2 )/2, 0), then ) Woo )<*foo = 0 • 111 mat case > algorithms 1 and 4 are 

f 

5 the same. This means that the term — a Woo )<* f oo is used to reduce the 

180 f f 2 -fi fi 

estimation error of the algorithm 1 resulting from non-symmetry of the phase response 
♦(f). 

For the real measured phase response, <|>(f) is only given at the sampling points 
of the frequency f. That is, (|>(f) is a discrete sequence rather than a continuous function 

1 0 in the interval [t, , f 2 ) . Furthermore, tff , ) and *(f 2 ) may be not equal to 1 80 and -1 80, 
respectively, although it is true that +(' 1 - Af ) < 0 , +(f j ) > 0 , tff 2 ) < 0 and tff 2 + Af ) > 0 , 
where Af is the sampling interval in the frequency domain. In order to increase the 
estimation accuracy, Af should be as small as possible and linear regression can be used 
to extrapolate 4>(f) from [f 1f f 2 )to [f n ewi» f new2)such that Wnewi) - 180 and 

1 5 Wnewi - 0) = -180 . For the extrapolated <K0, 

fo f new2 ~ W1 360 f 0 f 0 180 '0 T new2 ~ T new1 Uwt 



or simply 



*0 f new2 " T new1 360 f 0 T 0 



Notably, the above algorithms for estimating the propagation velocity are only 
20 based on a proper segment of the measured frequency response of the SH-SAW sensor. 
In at least this way, the algorithms are different than conventional methods. 
Conventional methods, may measure the time delay of the sensor or time delays 
between several fingers of IDT by some time domain device. Alternatively, 
conventional methods may be based on an entire phase frequency response, an 
25 amplitude frequency response, and an inverse Fourier transform. In contrast, the 

techniques described herein are based only on a segment of the measured frequency 



i 
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response of the SH-SAW sensor, i.e., the proper segment, which can be identified as 

described herein. • 

Algorithm 1 was implemented in the software platform Lab View for on-line 
velocity measurement. LabView is a commercially available software program that can 
5 be obtained from National Instruments, Inc., USA. Using a SAW sensor as a detector 
is generally based on the comparison of propagation characteristics of the sensor 

without and with a particular surface perturbation. Such a comparison can be carried 

> 

out dynamically or statically. 

For the dynamic comparison, the propagation characteristic of the sensor is 

10 measured as a time series. For the static comparison, the sensor should consist of two 
channels, the reference channel and detection channel (also called the active channel). 
For example, the reference channel may correspond to the reference portion between 
IDTs 15 and 17 (FIG. 1), whereas the detection channel may correspond to the active 
portion between IDTs 14 and 16. The surface perturbation takes place in the detection 

1 5 channel. In that case, the propagation characteristics of two channels of the sensor are 

measured. The measurement of propagation characteristic of SH-SAW sensor is the 
base for using the sensor as a detector. 

Resonance frequency f rC sonancc and propagation velocity v are two of the most 
direct propagation characteristics, and Av/v are more sensitive than Af/f to the surface 

20 perturbation. Therefore, the measurement of the propagation velocity may provide an 
improved method for efficiently using SH-SAW sensor as a detector. In contrast, 
conventional techniques typically measure the log-amplitude response A(f) and phase 
response <|>(f) of SH-SAW sensor by a sensor analyzer and read out the phase at the 
running frequency from the measured phase frequency response. Developing routines 

* 

25 for the direct measurement of the propagation velocity of the SH-SAW sensor becomes 
very important in the practical application of SH-SAW sensor. The algorithms 
presented above provide a software-based approach to estimate the propagation 
velocity from a proper segment of the measured phase frequency response, e.g., which 
can be measured by an 8753ET network analyzer commercially available from Agilent 

30 Technologies, Inc., USA. Other similar network analyzers could also be used. 

LabView was chosen as any easy way to both communicate with the sensor 
analyzer and to do the calculation based on the model. An exemplary LabView screen 
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associated with execution of the techniques described herein is shown in FIG. 5. The 
inputs are in part defined by the spectral resolution of the sensor analyzer and include 
the minimum and maximum frequency (equal to a total of 1 600 data points on the 
sensor analyzer) and the running frequency. 

* 

5 The outputs from Lab View allow various interpretations of the results to be 

made. For the majority of the plots shown on the Lab View screen shown there are two 
outputs, one representing the measurement sensor and the other representing the 
reference sensor. In Lab View, plots of phase and velocity can be used to determine if a 
sensor is working reliably. If the phase becomes noisy or if the amplitude decreases 
10 more than about 20db then the sensor is deemed bad and discarded. 

The phase histogram plot can be used as another measure of sensor quality and 
relates to non-linearity in the phase plot. The velocity plot may provide the result of 
the estimation according to the techniques described herein. 

Plots of magnitude and phase are taken at the running frequency and plotted 
1 5 against time. Lastly, temperature and analog values are outputs that are plotted against 
time. In order to calculate velocity from a proper segment of phase frequency response 
an internal sub- VI is used within Lab View. The first step is to calculate all the +180 
and -180 degree phase points. In other words, to determine the phase inflection 
frequencies, the technique may sample a plurality of phase responses at frequencies 
20 proximate to the running frequency and estimate the phase inflection frequencies as a 
function of the plurality of phase responses. 

The difference between the velocity expression in continuous frequency and its 
discrete realization was discussed above. However, this first step may also be applied 
during the implementation of Algorithm 1 in LabView. For example, if the sensor 
25 analyzer has a limited resolution in frequency, the data about a -1 80 to 1 80 degree 
phase transition may be similar to that below: 

-175.97, -176.866, -177.761, -178.657, -179.552, 179.5522, 178.6567, 
177.7612, 176.8657, 175.9701, 

An extrapolated linear fit may be used to find the exact -180 and 180 degree 
30 phase points. An illustration of exemplary + 1 80 and -1 80 degree phase points is 

shown in FIG. 6. FIG. 7 illustrates use of a bilaterally extrapolated linear fit to find the 
exact -180 and 180 degree phase. As can be appreciated from FIG. 7 it is possible to 
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* 

calculate the values of ft, fi and the phase at the running frequency <|>(fo). However, 
because the format is different in LabView it may be necessary to use different 
nomenclature from the model derived above. For example, in LabView: 

Tau_Fr = (((Fl/(F2-Fl))*360)-Pha_Fr+180)/F_ran*1000000*360) (with 
5 Fjun converted to hertz), 

and then, 

V_Fr=(L*0.001)/Tau_Fr (with L converted from millimeters to meters). 
The final output of velocity (V_Fr) is plotted in LabView and corresponds to v(f 0 ) , 

* 

discussed above. 

10 Using a separate application of LabView, it is possible; to take stored phase files 

and reconstruct velocity data. The phase data can then be transformed using exactly the 
same model as above. In this way, it is possible to see the effect of changing the 
running frequency relative to the velocity obtained. LabView allows two running 
frequencies (at 0.1MHz difference) to be evaluated simultaneously. Advantageously, 

15 the phase inflection frequencies define edges of a monotonically changing subset of a 
graph of phase versus frequency of the surface acoustic wave sensor. Thus, the 
techniques described herein may allow not only for the identification of a material via 
the sensor, but also an indication of the concentration of the material. 

20 Example 1 

The Love mode SH-S A W sensor used in the experiment was a LiTa03 device 
operating at 103 MH 2 provided by Sandia National Laboratories, USA. A low-walled 
flow cell was placed over the sensor and filled with Phosphate Buffered Saline (PBS) 
buffer solution at pH 7.5. This liquid container was connected to a syringe pump 

■ 

25 system to allow a slow flow of buffer. During the experiment, multiple aliquots of 250 
microliters of 0.05 mg/ml Bovine Serum Albumin %(BSA) protein was injected into 
the cell at designated times. An 8753ET network analyzer from Agilent Technologies, 
Inc., USA, measured the log-amplitude and phase response of the sensor approximately 
every twenty seconds. Based on the resulting phase responses, the curve of phase at the 

30 operation frequency versus time was immediately obtained. Algorithm 1 was then used 
to calculate the propagation velocity versus time at the operating frequency. FIGS. 8 
and 9 respectively show the phase response and propagation velocity plotted over the 
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time of the experiment. The valleys of the curve of velocity versus time precisely 

■ 

correspond to moments when BSA was injected into the very slowly flowing buffer 
stream. Compared with the curve of phase versus time, the calculation of propagation 
velocity provides much more information about the mass loading on the surface of a 
Love mode SH-SAW sensor, resulting in an increase in sensitivity to mass loading and 
viscosity. 



Example 2 

In this example, three separate concentrations of bacteria being detected were 
10 used. The first bacterial concentration was 10 3 cfu/ml (colony forming unit per 
milliliter), and the results are illustrated in FIGS. 10 and 1 1. The second bacterial 
concentration was 10 5 cfu/ml and the results are illustrated in FIGS. 12 and 13. The 
third bacterial concentration was 10 7 cfu/ml and the results are illustrated in FIGS. 14 
and 15. 

AV 

15 In particular, the graphs of FIGS. 10, 12 and 14 show— versus time at the 

specified concentrations, and the graphs of FIGs. 11,13 and 15 show^ versus time 

respectively at the specified concentrations. The usual way to use an SH-SAW sensor 
for mass-loading measurement is that the sensor contains two channels: reference 
channel and detection channel (the detection channel also being referred to herein as 

20 the "active channel" or "active portion" of the sensor). In the reference channel, the 
surface of the sensor is largely unperturbed, while it is perturbed in the detection 
channel. By measuring and comparing the propagation characteristics of the two 
channels, the perturbation resulting from the mass-loading in the detection channel can 
be detected. If the indicator function, which can be used to describe the difference of 

25 the propagation characteristic of the two channels, is monotonic with respect to the 
mass of the loading material, then the existence of loading mass can be detected and 

Av 

quantified. It has been pointed out above that using — as indicator function for 

detection is better than using — . This example demonstrates the correctness of this 
conclusion from a real experiment. In this experiment, the surface of the detection 
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channel was provided with bacterial loading at specific concentration. The 

3 5 7 

concentrations of the bacteria used in the experiment are 10 , 10 , and 10 cfu/mL, 
respectively. The delay from injecting bacteria into the detection channel until back- 
washing was around 20 minutes. The curves of — and versus time for each 

v 0 <(>o 

5 bacterial concentration are shown in FIGS. 10-15, where the calculations of the 

propagation velocities of the two channels are based on the above algorithm. It should 
be pointed out that the experiment in this example includes no* antibody. Bacteria were 
bound directly to the surface chemistry which was a saccharin immobilization layer. 

The — and — at the steady state are extracted from the above curves for each 
v 0 *o 

■ 

4 

10 bacterial concentration. The relationships of — versus log-concentration and 

v o 

^ versus log-concentration are respectively shown in FIGS. 16 and 17. 

• Although the experimental results are just for three different bacterial 
concentrations, clearly, ^ is not monotonic with respect to the bacterial concentration. 

It also demonstrates the correctness of the conclusion that ^ is not a proper indicator 

<Po 

1 5 for quantitatively detecting loading mass with large dynamic range. 



Example 3 

The Love mode SH-SAW sensor used in the experiment was a LiTa0 3 device 
operating at 103 MH Z provided by Sandia National Laboratories, USA. During the 
20 experiment, the surface of the sensor was: 

1 . In air for approximately 20 min 

2. In buffer for approximately 30 mih 

3. Subjected to bacteria in the buffer for approximately 18 min 

4. Briefly washed and left in buffer for approximately 20 min 
25 5. Subjected to De Ionized (DI) water for approximately 30 min. 

Bacteria were bound to the antibody attached to the surface. FIG. 18 is a graph 
of the propagation velocity at the operating frequency versus time and FIG. 19 is a 
graph of phase at the operating frequency versus time. The propagation velocities are 
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calculated using algorithm 1 presented above. In FIGS. 18 and 19, the dashed vertical 
lines define the time durations of various surface conditions of the sensor as detailed 
above. 

By observing the graph of velocity versus time in FIG. 18, the process and 

5 results of the experiment can be understood much more easily. The propagation of the 

surface acoustic wave of the sensor has higher velocity when its surface was in air as 

shown at 181. Upon injecting buffer, the velocity decreased and tended to a steady 

■ 

state as shown at 182. Upon adding a certain amount of bacteria, the velocity 
decreased again as shown at 1 83. The brief washing did not remove the bacteria, which 

10 were adsorbed by the antibody attached to the surface of the sensor. Therefore, the 
velocity during washing as that in the previous state as shown at 184. Finally, the 
velocity of the surface acoustic wave sensor increased because the surface condition 
was changed from the adsorbed bacteria to water as shown at 185. By stark contrast, it 
is difficult to understand and explain the experimental process and results based on the 

1 5 curve of phase versus time shown in FIG. 1 9. 

Example 4 

The Love mode SH-S AW sensor used in the experiment was a LiTa03 device 
operating at 103 MH Z provided by Sandia National Laboratories, USA. A low-walled 

20 flow cell was placed over the sensor and filled with Phosphate Buffered Saline (PBS) 
buffer solution at pH of approximately 7.5. This liquid container was connected to a 
syringe pump system to allow a slow flow of buffer. During the experiment, multiple 
aliquots of 250 microliters of Bovine Serum Albumin (BSA) protein at various 
concentrations were injected into the cell at designated times. The injection of proteins 

25 conformed to Table 1, below. In Table 1, BSA refers to Bovine Serum Albumin, and 
IgG refers to Immunoglobulin. 
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TABLE1 



Injection 

moment 

(sampling 

points in 

time 

domain) 


Injected 
Proteins 


Concentration 
(rag/ml) 


Age of Solution • 

t 

• 

• 


Locations of 
jumping 
down point 
of velocity 


Av/v 


70 


BSA 


0.05 


1 day old 


71 


-0.0193 


96 


IgG 


0.05 


2 weeks old 






122 


BSA 


0.05 


1 day old 


126 


-0.0196 


159 


Buffer 




• 






184 

» 


BSA 


0.05 


2 weeks old 






737 


BSA 


0.05 


1 day old 


741 


-0.0197 


783 


BSA 


0.025 


1 day old 


787 


-0.0162 


816 

* 


BSA 


0.0125 


1 day old 


820 


-0.01 10 


854 


BSA 


0.0333 


1 day old 


858 


-0.0169 


888 


BSA 


0.04 


1 day old 


892 


-0.01 86 


916 


BSA 


0.05 


1 day old 


920 


-0.0201 


948 


BSA 


0.00625 


1 day old 


952 


-0.0056 


994 


BSA 


0.003125 


1 day old 


997 


-0.0011 


1023 


BSA 


0.05 


1 day old 


1027 


-0.0201 


1058 


IgG 


0.05 


1 day old 


1061 


-0.0014 


1084 


BSA 


0.05 


1 day old 


1087 


-0.0209 

■ 



As described in Example 1, an 8753ET network analyzer from Agilent 
Technologies, Inc., USA, measured the log-amplitude and phase frequency response of 
5 the sensor approximately every 20 seconds. Based on proper segments of these phase 
frequency responses for the active channel, and algorithm 1 proposed above, both the 
curves of phase and propagation velocity at the running frequency versus time fare 
were obtained. FIG. 20 is the resultant graph of phase versus time, whereas FIG. 21 is 
the resultant graph of propagation velocity versus time. 
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The details of the curve of propagation velocity (m/s) versus time (sampling 
point in time domain) for BSA injections are shown in the Figures 22-24. The locations 
of the valleys of the curve are very close to the injection moments listed in Table 1. In 
FIG. 23 the labels correspond to the following: 
5 (1) (737) BSA .05 mg/ml; 

(2) (783)BSA.025mg/ml; 

(3) (816) BSA .0125 mg/ml; 

(4) (854) BSA .0333 mg/ml; 

(5) (888) BSA .04 mg/ml; 
10 (6) (916) BSA .05 mg/ml; 

(7) (948) BSA 0.00625 mg/ml; 

(8) (994) BSA 0.003125 mg/ml; 

(9) (1023) BSA .05 mg/ml; 

In order to estimate Av/v more accurately, some segments of the curve of 

1 5 velocity around the injection moments of BSA are used. An example is shown in FIG. 
26. In this case, the estimated injection moment is the 952 nd sampling point in time 
domain. From this segment illustrated in FIG. 26, the edges of the valley are 
determined and the velocities v 0 and Vi are estimated. Then, the Av/v=(V|-v 0 )/v 0 for 
this injection is calculated, in this case, Av/v=-0.0056. FIG. 25 shows the 

20 corresponding segment of the curve of <j>. It is difficult to determine A<|>/<J> for this 
injection from the graph of FIG. 25. 

More examples are presented in FIGS. 27-30. In particular, FIGS. 28 and 30 
illustrate segments of the curve of v around the estimated injection moments and FIGS. 
27 and 29 illustrate the corresponding segments of the curve of <J>. The calculated 

25 (Av/v)'s for FIGS. 28 and 30 are -0.0201 and -0.0197, respectively. It is still difficult 
to determine A<M> from FIGS. 27 and 29. 

The more accurate estimations of Av/v about the estimated injection moments 
of BSA are presented in Table 1. The resulting Av/v versus the known concentration of 
the injected BSA at the same injection moment, which is listed in Table 1, is shown in 

30 FIG. 31 as a circle. Corresponding to the concentration 0.05 mg/ml, there are seven 
overlapping (Av/v)'s.° 
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One can examine the circle corresponding to the largest Av/v. The curve of 
FIG. 31 results from connecting the circles. The curve of FIG. 31 is monotonic. 

» 

Moreover, the monotonic nature of the curve of FIG. 31 demonstrates that using Av/v, 
which can be determined by the curve of velocity and calibration, the sensor may 
5 further detect the concentration of the injected BSA if the concentration of the injected 

BSA is in the interval [0.003125, 0.05] mg/ml. 

The measurement of propagation characteristic of SH-SAW sensor is the base 
for using the sensor as a detector. The resonance frequency f rC sonancc and propagation 
velocity v are two direct propagation characteristics and Av/v is more sensitive than 

1 0 Af/f to the surface perturbation. Therefore, the measurement of the propagation 
velocity is a preferred index for efficiently using SH-SAW sensor as a detector. 

In contrast, the conventional approach is to measure the log-amplitude 
frequency response A(f) and phase frequency response (j)(f) of SH-SAW sensor by a 
sensor analyzer, and then read out the phase at the running frequency. The serious 

1 5 limitation of A<|> or AM for detection was discussed above. By developing hardware 
and/or software for the direct measurement of the propagation velocity of the SH-SAW 
sensor, a practical application of a SH-SAW sensor can be achieved. Based on the 

h 

system model of SH-SAW sensor with conventional triple transit echo (TTE) 
techniques, several algorithms have been presented to estimate the propagation velocity 
20 from a proper segment of the phase frequency response measured by a sensor analyzer. 

■ 

The first algorithm (algorithm 1) has also been implemented in the software platform 
Lab View for on-line velocity measurement. 

Compared with the conventional phase-based or phase shift- based method, the 
techniques described herein will not lose any information carried by phase. Also, 
25 because of overcoming the phase ambiguity, the techniques described herein can 

extract much more information about the surface condition of SH-SAW sensor relative 

■ 

to conventional techniques. Furthermore, the examples showed the monotonic nature 
of the curve of Av/v versus concentration of the liquid mass loading in a large dynamic 
range of the concentration. This may therefore be the basis for using SH-SAW sensor 
30 in practice to quantitatively detect mass loading. In contrast, for the curve of AM 
versus concentration of the liquid mass loading, monotonic curves only hold true for 
small concentrations. 
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Various embodiments of the invention have been described. In particular, 
techniques for estimating the propagation velocity through a surface acoustic wave 
sensor have been described. The techniques may be implemented in hardware, 
software, firmware, or the like. Example hardware implementations include 
5 implementations within a general purpose microprocessor, an application specific 

integrated circuit (ASIC), a field programmable gate array (FPGA), specifically 
designed hardware components, or any combination thereof. In addition, one or more 
of the techniques described herein may be partially or wholly executed in software. In 
that case, a computer-readable medium may store or otherwise comprise computer- 

10 readable instructions, i.e., program code, that can be executed. by a processor or to carry 
out one of more of the techniques described above. The techniques can be used for 
both constant time delay (in non-dispersive case) and frequency dependant time delay 
(in dispersive case). 

For example, the computer-readable medium may comprise random access 

1 5 memory (RAM), read-only memory (ROM), non- volatile random access memory 

(NVRAM), electrically erasable programmable read-only memory (EEPROM), flash 
memory, or the like. These and other embodiments are within the scope of the 
following claims. 

Suitable methods for coating the devices of the present invention include 
20 Applicants' Copending Application U.S. Serial No. 10/607,698, filed June 27, 2003. 

The present invention may be utilized in combination with various materials, 
methods, systems, apparatus, etc. as described in various U.S. patent applications 
identified below, all of which are incorporated by reference in their respective 
entireties. They include: U.S. Patent Application Serial Nos. 60/533,162, filed on 
25 December 30, 2003; 60/533,178, filed on December 30, 2003; 10/896,392, filed July 
22, 2004; 10/713,174, filed November 14, 2003; 10/987,522, filed November 12, 2004; 
10/714,053, filed November 14, 2003; 10/987,075, filed November 12, 2004; 
60/533,171, filed December 30, 2003; 10/960,491, filed October 7, 2004; 60/533,176, 

filed December 30, 2003; 60/533,169, filed December 30, 2003; , titled 

30 "Method of Enhancing Signal Detection of Cell- Wall Components of Cells", filed on 

even date herewith (Attorney Docket No. 59467US002); , titled "Soluble 

Polymers as Amine Capture Agents and Methods", filed on even date herewith 
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(Attorney Docket No. 59995US002); , titled "Multifunctional Amine 

Capture Agents", filed on even date herewith (Attorney Docket No. 59996US002); 

PCT Application No. , titled "Surface Acoustic Wave Sensor Assemblies", 

filed on even date herewith (Attorney Docket No. 58928WO003); PCT Application No. 

5 , titled "Acousto-Mechanical Detection Systems and Methods of Use", 

filed on even date herewith (Attorney Docket No. 59468 WO003); PCT Application No. 

, titled "Detection Cartridges, Modules, Systems and Methods", filed on 

even date herewith (Attorney Docket No. 60342 WO003); and PCT Application No. 
9 titled "Acoustic Sensors and Methods", filed on even date herewith 

10 ' (Attorney Docket No. 60209 WO003). 

The complete disclosures of the patents, patent applications, patent documents, 
and publications cited herein are incorporated by reference in their entirety as if each 
were individually incorporated. Various modifications and alterations to this invention 
will become apparent to those skilled in the art without departing from the scope and 

1 5 spirit of this invention. It should be understood that this invention is not intended to be 
unduly limited by the illustrative embodiments set forth herein and that such 
embodiments are presented by way of example only, with the scope of the invention 
intended to be limited only by the claims. 
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CLAIMS: 

1 . A method comprising: 

identifying a segment of phase frequency response of a surface acoustic wave 

5 sensor; and 

estimating a time delay associated with wave propagation through the surface 
acoustic wave sensor based on the identified segment of phase frequency response. 

2. The method of claim 1, wherein identifying the segment of phase frequency 
1 0 response comprises determining first and second phase inflection frequencies 

proximate to a running frequency associated with the surface acoustic wave sensor. 

3. The method of claim 2, wherein determining phase inflection frequencies 
comprises: 

1 5 sampling a plurality of phase responses at frequencies proximate to the running 

frequency and initially estimating phase inflection frequencies as a function of the 
plurality of phase responses at frequencies proximate to the running frequency; 

sampling a plurality of phase responses at frequencies proximate to the initially 
estimated phase inflection frequencies; and 

20 more accurately estimating the phase inflection frequencies as a function of the 

plurality of phase responses at frequencies proximate to the initially estimated phase 
inflection frequencies. 

4. The method of claim 2, wherein the first and second phase inflection 

25 frequencies define edges of a monotonically changing subset of a graph of phase versus 
frequency of the surface acoustic wave sensor. 

5. The method of claim 2, further comprising estimating the time delay according 
to approximately the following equation: 

30 LM t M 

/./■-/, 360 /„ /, 
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where f(f 0 ) is the time delay at frequency /„ , /„ is the running frequency, f, is the 
first phase inflection frequency, f 2 is the second phase inflection frequency, and 
4>(f 0 ) is a measured phase response of the surface acoustic wave sensor at frequency 

fo- 

5 

6. The method of claim 1 , further comprising estimating the time delay according 

i 

to approximately the following equation: 

where r (/ 0 ) is the time delay, f 0 is a running frequency, <K/ 0 ) is a measured phase 
10 response of the surface acoustic wave sensor, /♦ is any frequency between a first phase 
inflection frequency and a second phase inflection frequency, ^(/.) is a measured 
phase frequency response at the frequency /• , and <*(/•) is a first order of derivative of 
the measured phase response at the frequency f. . 

15 7. The method of claim 1 , further comprising estimating the time delay according 
to approximately the following equation: 

* 

where f (/ 0 ) is the time delay, and ^(/ 0 ) is a first order of derivative of a measured 

< 

phase response at a frequency f 0 . 



20 



8. The method of claim 1 , further comprising estimating the time delay according 
to approximately the following equation: 

where f(/ 0 ) is the time delay, / 0 is an operating frequency, f y is a first phase 
25 inflection frequency, f 2 is a second phase inflection frequency, and ^(/ 0 ) is a measured 

h 

phase response of the surface acoustic wave sensor, integral \$( f 0 o)dfoo * s equal to 
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integral \<f>(f )df , where <p( f) is a measured phase response at frequency / and 
/ varies from / to f 2 . 



9. The method of claim 1, further comprising estimating a propagation velocity of 

* 

5 a surface acoustic wave through the surface acoustic wave sensor from an estimated 
time delay according to the following equation: 
L 

v(f) = , where v(f)is an estimated propagation velocity of the surface acoustic 

iff) 

wave at frequency / , f (f ) is the estimated time delay at the frequency / , and I is a 
distance between centers of an input inter-digitized transducer (IDT) and an output 
10 IDT. 



1 0. The method of claim 1 , further comprising: 

bringing a fluid into contact with a surface of the surface acoustic wave sensor; 

and 

1 5 identifying a concentration of a material in the fluid as a function of an 

estimated propagation velocity, the estimated propagation velocity being estimated 
based on the estimated time delay. 

11. The method of claim 1, wherein the surface acoustic wave sensor comprises a 
20 Love mode shear-horizontal surface acoustic wave sensor. 

12. A computer-readable medium comprising instructions that when executed in a 
processor: 

identify a segment of phase frequency response of a surface acoustic wave 
25 sensor; and 

estimate a time delay associated with wave propagation through the surface 
acoustic wave sensor based on the identified frequency response. 

13. The computer-readable medium of claim 12, further comprising instructions 
30 that when executed identify the segment of phase frequency response by determining 
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first and second phase inflection frequencies proximate to a running frequency 
associated with the surface acoustic wave sensor. 



14. The computer-readable medium of claim 13, further comprising instructions 
5 that when executed determine phase inflection frequencies for a discrete phase 

frequency response by: 

sampling a plurality of phase responses at frequencies proximate to the running 
frequency and initially estimating phase inflection frequencies as a function of the 
plurality of phase responses at frequencies proximate to the running frequency; 
10 sampling a plurality of phase responses at frequencies proximate to the initially 

estimated phase inflection frequencies; and 

more accurately estimating phase inflection frequencies as a function of the 
plurality of phase responses at frequencies proximate to the initially estimated phase 
inflection frequencies. 

15 

1 5. The computer-readable medium of claim 1 3, wherein the first and second phase 
inflection frequencies define edges of a monotonically changing subset of a graph of 

■ 

phase versus frequency of the surface acoustic wave sensor. 

20 1 6. The computer-readable medium of claim 1 3, further comprising instructions 
that when executed estimate the time delay according to approximately the following 

i' 

i 

equation: 

Vj-A_! L«1 ( oj 

hh~f\ 360 f 0 U 
where r (f 0 ) is the time delay at frequency / 0 , / 0 is the running frequency, f x is the 

25 first phase inflection frequency, f 2 is the second phase inflection frequency, and 

0(/ o ) is a measured phase response of the surface acoustic wave sensor at the running 
frequency f 0 . 



17. The computer-readable medium of claim 12, further comprising instructions 
30 that when executed estimate the time delay according to approximately the following 
equation: 



ii 
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where r(/ 0 ) is the time delay, f 0 is a running frequency, ^(/b) is a measured phase 
response of the surface acoustic wave sensor, f> is any frequency between a first phase 
inflection frequency and a second phase inflection frequency, ^(/*) is a measured 
5 phase frequency response at the frequency /. , and <p(f. ) is a first order of derivative of 
the measured phase frequency response at the frequency f> . 

18. The computer-readable medium of claim 12, further comprising instructions 
that when executed estimate the time delay according to approximately the following 
10 equation: 

where f (/ 0 ) is the time delay, and $(f 0 ) is a first order of derivative of a measured 

* 

phase response at a frequency f 0 . 

15 19. The computer-readable medium of claim 1 2, further comprising instructions 
that when executed estimate the time delay according to approximately the following 
equation: 

fofi-f) 360 Jo Jo 180 / 0 / 2 -/, j x 

where f (f 0 ) is the time delay, /„ is an operating frequency, f x is a first phase 
20 inflection frequency, f 2 is a second phase inflection frequency, (f>(f 0 ) is a measured 

h 

phase response of the surface acoustic wave sensor, integral \<f>( f w )df m is equal to 
h 

integral \$(f )df , where <f>(f) is a measured phase response at a frequency / and 
/. 

/ varies from /J to f 2 . 



II 
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20. . The computer-readable medium of claim 12, further comprising instructions 

■ 

that when executed estimate a propagation velocity of the surface acoustic wave from 
the estimated time delay according to the following equation: . 

v(f ) = - , where v(f)\s the estimated propagation velocity of the surface acoustic 
i(f) 

5 wave at frequency / , i (f) is the estimated time delay at frequency / , and L is a 

distance between centers of an input inter-digitized transducer . EDT and an output IDT. 

2 1 . The computer-readable medium of claim 1 2, further comprising instructions 
that when executed identify a concentration of a material in a fluid as a function of an 

1 0 estimated propagation velocity, the estimated propagation velocity being estimated 
based on the estimated time delay. 

22. The computer-readable medium of claim 12, wherein the surface acoustic wave 
sensor comprises a Love mode shear-horizontal surface acoustic wave sensor. 

15 

23. A system comprising: 

a surface acoustic wave sensor; 

a sensor analyzer to receive output of the surface acoustic wave sensor; and 
a processor to receive input from the sensor analyzer, identify a segment of 
20 phase frequency response of a surface acoustic wave sensor, estimate a time delay 

associated with wave propagation through the surface acoustic wave sensor based on 

the identified segment of phase frequency response. 

24. The system of claim 23, wherein the processor identifies the segment of phase 
25 frequency response by determining first and second phase inflection frequencies 

proximate to a running frequency associated with the surface acoustic wave sensor. 



25. The system of claim 24, wherein the processor determines the phase inflection 
frequencies by: 
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sampling a plurality of phase responses at frequencies proximate to the running 
frequency and initially estimating the phase inflection frequencies as a function of the 

■ 

plurality of phase responses at frequencies proximate to the running frequency; 

sampling a plurality of phase responses at frequencies proximate to the initially 
estimated phase inflection frequencies; and 

more accurately estimating the phase inflection frequencies as a function of the 
plurality of phase responses at frequencies proximate to the initially estimated phase 
inflection frequencies. 



10 26. The system of claim 24, wherein the first and second phase inflection 

frequencies define edges of a monotonically changing subset of a graph of phase versus 
frequency of the surface acoustic wave sensor. 



27. The system of claim 24, wherein the processor estimates the time delay 
1 5 according to approximately the following equation: 

f(/o) = A_J L£^) + H 

° foA-fi 360 f 0 f 0 
where f(fo) is ^ e time delay at frequency f 0 , f 0 is the running frequency, f x is the 
first phase inflection frequency, f 2 is the second phase inflection frequency, and 
$(fo) is a measured phase response of the surface acoustic wave sensor at the running 

* 

20 frequency / 0 . 



28. The system of claim 23, wherein the processor estimates the time delay 
according to approximately the following equation: 

,i 
i 

25 where r(/ 0 ) is the time delay, f 0 is a running frequency, <fi{f 0 ) is a measured phase 

response of the surface acoustic wave sensor, f, is any frequency between a first phase 
inflection frequency and a second phase inflection frequency, is a measured 
phase frequency response at the running frequency /. , and <f>{fj) is a first order of 
derivative of the measured phase frequency response at the running frequency f. . 
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29. The system of claim 23, wherein the processor estimates the time delay 
according to approximately the following equation: 

where f(/ 0 ) is the time delay, and ^(/ 0 ) is a first order of derivative of a measured 
phase response at frequency f 0 . 



30. The system of claim 23, wherein the processor estimates the time delay 
according to approximately the following equation: 

where r(/ 0 ) is the time delay, f Q is an operating frequency, f x is a first phase 

i 

« 

inflection frequency, f 2 is a second phase inflection frequency, and <f>{f Q ) is a measured 

* • 

phase response of the surface acoustic wave sensor, integral J <j>( f Q0 )df m is equal to 
h 

integral \<p( f )df , where <j>(f) is a measured phase response at frequency / ,and 

/. 

15 / varies from f } to / 2 . 



3 1 . The system of claim 23, wherein the processor estimates a propagation velocity 
of the surface acoustic wave based on the estimated time delay according to the 
following equation: 

20 Hf) = —— > where v( f)\s an estimated propagation velocity of the surface acoustic 

*(f) 

wave at a frequency / , r(f) is the estimated time delay at the frequency / , and L is 

■ 

a distance between centers of an input inter-digitized transducer DDT and an output 
IDT. 



25 32. The system of claim 23, wherein the processor estimates a propagation velocity 
based on the estimated time delay. 
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33. The system of claim 32, wherein the processor identifies a concentration of a 
material in a fluid as a function of the estimated propagation velocity. 

5 34. The system of claim 23, wherein the surface acoustic wave sensor comprises a 
Love mode shear-horizontal surface acoustic wave sensor. 
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